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Abstract. We consider a new type of defect in the scope of linear elasticity theory, 

using geometrical methods. This defect is produced by a spherically symmetric disloca- 
tion, or ball dislocation. We derive the induced metric as well as the affine connections 
and curvature tensors. Since the induced metric is discontinuous, one can expect am- 
biguity coming from these quantities, due to products between delta functions or its 
derivatives, plaguing a description of ball dislocations based on the Geometric Theory 
of Defects. However, exactly as in the previous case of cylindric defect, one can obtain 
some well-defined physical predictions of the induced geometry. In particular, we ex- 
plore some properties of test particle trajectories around the defect and show that these 
trajectories are curved but can not be circular orbits. 
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^ ■ 1 Introduction 

co 

^ The topological defects attract great interest due to the applications to condensed matter physics 

(see, e.g., [Tj [2] for an introduction and recent review). Another elegant approach in this area is 
the geometric theory of defects [31 H] (see also [5] for the introduction), which is formulated in 
terms of the notions originally developed in the theories of gravity. In this framework, we can cite 
basically two kinds of defects, described in the view of Riemann-Cartan geometry: disclinations 
and dislocations. This means that the curvature and torsion tensors, respectively, are interpreted as 
surface densities of Frank and Burgers vectors and thus linked to the nonlinear, generally inelastic 
deformations of a solid. Recently, the qualitatively new kind of geometric defect has been described 
in [6] (see also [7]), corresponding to a tube dislocation (with cylindrical symmetry). 

Nevertheless, one can find others well known types of dislocations. In Elasticity Theory, the 
dislocations and their physical effects are a matter of great concern, specially the screw dislocation 
[8]. In this paper, we consider a ball dislocation (or sphere dislocation), which is the same type of 
defect studied in [6] , translated to spherical symmetry. Is is remarkable that such a problem was not 
investigated yet. For our purposes, the approach will be limited to linear elasticity theory (using 
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Riemann-Cartan geometry as a tool). The ball dislocation can be understood by the illustration 
in Figure 1. 




Figure 1: Ball dislocation produced by cutting a spherical sector, and then indentifying the surfaces 
r\ and T2 with p. 



It is possible to treat this defect as a point defect in the limiting case where p is much smaller 
than the dimensions of the considered physical system. Of course, the ball dislocation provides 
a more general picture. This paper is organized in the following way. In Section 2, we treat 
mathematically the ball dislocation, find the solution of the variable which describes the defect, 
and with the help of Gravity Theory methods, we study the trajectories of test particles in Section 
3. Finally we proceed to our conclusions. 



2 Ball dislocation in linear elasticity theory 

Let us describe the ball dislocations in linear elasticity theory. Consider a homogeneous and 
isotropic elastic media as a three-dimensional Euclidean space M 3 with Cartesian coordinates x % , y l , 
where i = 1,2,3. The Euclidean metric is denoted by 5ij = diag(+ + +). The basic variable in 
the elasticity theory is the displacement vector of a point in the elastic media, u l {x), x € M 3 . In 
the absence of external forces, Newton's and Hooke's laws reduce to three second order partial 
differential equations which describe the equilibrium state of elastic media (see, e.g., [9]), 

(1 - 2a)A Ui + didju j = 0. (1) 

Here A is the Laplace operator and the dimensionless Poisson ratio a (—1 < a < 1/2) is defined 
as 



a 



A 
2(X + p) 



The quantities A and p, are called the Lame coefficients, which characterize the elastic properties 
of media. 

Raising and lowering of Latin indices can be done by using the Euclidean metric, Sij, and its 
inverse, 5 tJ . Eq. ([1]) together with the corresponding boundary conditions enable one to establish 
the solution for the field u l {x) in a unique way. 

Let us pose the problem for the ball dislocation shown in Figure 1. This dislocation can be 
produced as follows. We cut out the thick spherical sector of media located between two concentric 
spherical surfaces of radii r\ and T2 (t\ < r-z), move symmetrically both cutting surfaces one to the 
other and finally glue them. Due to spherical symmetry of the problem, in the equilibrium state 
the gluing surface is also spherical, of radius p which will be calculated below. 

Within the procedure described above and shown in Figue 1, we observe the negative ball 
dislocation because part of the media was removed. This corresponds to the case of r\ < r%. 
However, the procedure can be applied in the opposite way by addition of extra media to M 3 . In 
this case, we meet a positive ball dislocation and the inequality has an opposite sign, T\ > r%- 

Let us calculate the radius of the equilibrium configuration, p. This problem is naturally formu- 
lated and solved in spherical coordinates, r, 8, <f>. Let us denote the displacement field components 
in these coordinates by u r , u s , vF. In our case, u 8 = u^ = due to the symmetry of the problem, 
so that the radial displacement field u r (r) can be simply denoted as u r (r) = u(r). 

The boundary conditions for the equilibrium ball dislocation are 



duh 



u\ „ = , u\ = , 

lr=0 lr=oo A j. 



du e 



* - < 2 > 

The first two conditions are purely geometrical, and the third one means the equality of normal 
elastic forces inside and outside the gluing surface in the equilibrium state. The subscripts "in" 
and "ex" denote the displacement vector field inside and outside the gluing surface, respectively. 

Let us note that our definition of the displacement vector field follows [5] , but differs slightly 
from the one used in many other references. In our notations, the point with coordinates y' 1 , after 
elastic deformation, moves to the point with coordinates x 1 : 

y^ X \y)=tf + U l (x). (3) 

The displacement vector field is the difference between new and old coordinates, u l {x) = x % — y l . 
Indeed, we are considering the components of the displacement vector field, u l (x), as functions of the 
final state coordinates of media points, x l , while in other references they are functions of the initial 
coordinates, y l . The two approaches are equivalent in the absence of dislocations because both 
sets of coordinates x l and y l cover the entire Euclidean space M 3 . On the contrary, if dislocation is 
present, the final state coordinates x l cover the whole M 3 while the initial state coordinates cover 
only part of the Euclidean space lying outside the thick sphere which was removed. For this reason 
the final state coordinates represent the most useful choice here. 

The elasticity equations ([1]) can be easily solved for the case of ball dislocation under considera- 
tion. Using the Christoffel symbols in order to evaluate the expressions for the differential operators 
(Laplacian and divergence), equations ([T]) reduce to the only one non-trivial equation 

d 2 u{r) 2 du(r) 2 
dr 2 r dr r 2 



One can remember that only the radial component differs from zero. The angular 9 and <j) compo- 
nents of equations (P) are identically satisfied. The general solution for ^ is given by 

U = aT ~ ^2 ' ( 5 ) 

which depends on the two arbitrary constants of integration a and /?. Due to the first two boundary 
conditions ([2]), the solutions inside and outside the gluing surface are 

u- m = ar , a > 0, 

P (6) 

^ex = 2 ' /5 > 0. 

The signs of the integration constants correspond to the negative ball dislocation. For positive ball 
dislocation, both integration constants have opposite signs. 

Using the solution ([5]) and the third boundary condition ([2]), one can determine the radius of 
the gluing surface, 

»=^i. (7) 

One can see that p is not the mean between r\ and 7*2, as it is for the cylindrical symmetry defect 
[6]. On the contrary, the gluing surface is located closer to the external radius 7*2. After simple 
algebra, the integration constants can be expressed in terms of p and the thickness of the sphere 

I = r 2 — 7*1 : 

21 p 2 l 



a = s~ P » = T- (8) 



It is straightforward also to get 



3/9-2/ , 3p + l 

7*1 = and ro = . 

3 3 

Observe that as r% is positive, we must have always p > 2Z/3. 

Finally, within the linear elasticity theory, eq. ([6]) with the integration constants (JSj) yields a 
complete solution for the ball dislocation in linear elasticity theory, valid for small relative displace- 
ments, when l/n <C 1 and ljr% <C 1. It is remarkable that the solution obtained in the framework 
of linear elasticity theory does not depend on the Poisson ratio of the media. In this sense, the ball 
dislocation is a purely geometric defect which does not feel the elastic properties. 

In order to use the geometric approach, we compute the geometric quantities of the manifold 
corresponding to the ball dislocation. From the geometric point of view, the elastic deformation 
is a diffeomorphism between the given domains in the Euclidean space. The original elastic media 
M 3 , before the dislocation is made, is described by Cartesian coordinates y l with the Euclidean 
metric Sij. An inverse diffeomorphic transformation x — >• y induces a nontrivial metric on M 3 , 
corresponding to the ball dislocation. In Cartesian coordinates, this metric has the form 

dy k dy l 
9 ^ x) = ^d^J 6kl - (9) 

We use curvilinear spherical coordinates for the ball dislocation and therefore it is useful to modify 
our notations. The indices in curvilinear coordinates in the Euclidean space IR 3 will be denoted 



by Greek letters x* 1 , jjl = 1,2,3. Then the "induced" metric for the ball dislocation in spherical 
coordinates is 



QyP QyV 

dx^ dx v 



aA x ) = —,^ 9 P a, (io) 



where g pu is the Euclidean metric written in spherical coordinates. We denote spherical coordinates 
of a point before the dislocation is made by {y, 9,<p}, where y without index stands for the radial 
coordinate and we take into account that the coordinates 9 and <j) do not change. Then the 
diffeomorphism is described by a single function relating old and new radial coordinates of a point, 
y = r — u(r), where 

u{r) = —r , r < p (11) 

dp 



pH 

3r 2 



u(r) = Hb> r >P- ( 12 ) 



It is easy to see that this function has a discontinuity -u ex t — u- m t\ r =p = I at the point of the cut. 
Therefore a special care must be taken in calculating the components of induced metric. 

It is useful to express u(r) in a way simultaneously valid in both domains, r < p and r > p. We 
have then 



u(r) = —rH - r) - -^ H(r - p) , (13) 

where H(r — p) is the Heaviside step function. As H'{r — p) = 5(r — p), one achieves 

= ^M = v , r ) _ ls , r _ ) (14) 

dr 

where 

21 . s 2p 2 l . 

By direct calculation of induced metric, by (|10p . one can write the corresponding line element 

as 

ds 2 = (1 - u'fdr 2 + (r - u) 2 (d9 2 + sin 2 6d0 2 ) . (15) 

It is clear that the above expression, besides discontinuous, contains also a (^-function thanks to v! . 
In order to avoid further conceptual consequences coming from a 5-function in the line element, 
we shall drop it, and adopt u'(r) —> v(r). In other words, let us consider the line element 

2l\ 2 

— i 1 

hit 



ds 2 nt - (1-— ) (dr 2 + r 2 d9 2 + r 2 sin 2 9d<j) 2 ) , 

ds 2 ext = (l - |J) 2 dr 2 + (l + 0) \r 2 d9 2 + r 2 sin 2 9dtf) . 

Notice that the interior space is conformally flat, with a constant scale factor, while the exterior 
metric is not so. Both metrics are flat (as follows by direct calculation of Riemann tensor), as they 



^^With a 5-function in the metric, the Burgers vector can not be defined properly as a surface integral [5]. 



should be (because they were obtained by coordinate transformations starting from the Euclidean 
metric). Nevertheless, the whole space is non-trivial since curvature is non-trivial exactly in the 
gluing surface. Next, we are going to investigate the consequences for trajectories of test particles 
around the defect. 

3 Trajectories of test particles around the defect 

What are the trajectories of test particles in a space with such a defect? Of course, the trajectories 
without defect would be straight lines, so we expect deviation from straight lines in the actual path. 
How these trajectories can be described? Are there any possible closed path around the defect? 
In order to answer these questions, let us consider the geodesic equations for both metrics, in the 
interior and in the exterior of gluing surface. The geodesic equations read 

— + r» pX uw x = o, (16) 

where U^ = dx^/dr = (1, r, 6, <j>) (c = 1 and dot means derivative with respect to r) and 

rM pA = -jQ^ {dxgap + d p g a \ - d a g pX ) ■ 

Let us remmember that if only dislocations are present, then only torsion (without Rieman- 
nian curvature) is found in the geometric approach - only the Burgers vector is non-trivial. But 
for practical purposes, one can treat the problem in the reverse way, considering only Rieman- 
nian curvature, because both approaches are equivalent. This equivalence is very well-known in 
telleparalelism (see, e.g., [10J). 

Inside the defect, the geodesies are straight lines and thus we shall consider only the exterior 
metric. By direct calculation, the geodesic equations, outside the gluing surface, can be written as 

6p 2 / , 2 <^l + 3rf)/- 22 - 2 \ 

4^2; _ c 3 

-r#-sin0cos^ 2 = O, (18) 



r(p 2 l + 3r 3 ) 

4o 2 Z-6r 3 ■ 2cos(9-- 
r(p z l + 3r d ) sinfl 

The denominator appearing on (JTTj) is always positive, because from condition p > 2Z/3, one gets 
3r 3 - 2p 2 l > 3r 3 - 3p 3 > 0. 

An interesting feature that we can understand is that, if r = 0, then the test particle should 
be necessarily at rest. This follows from equation ()17p . This means that if the test particle is 
moving, so its radial coordinate must be changing: there is no possible trajectory confined in a 
spherical surface. In other words, one can say that such a geometrical defect can not serve as an 
alternative description of gravitating objects (around which we know there are permitted circular 



orbits). However, other kinds of effects, in condensed matter physics, for example, can not be ruled 
out. 

A test particle can follow also a radial path, defined as any trajectory described by 9 = 4> = 0. 
To see that, let us consider 6 = cf> = 0, such that the radial equation (fT7|) reads 

6p 2 l 



r 2 



which can be solved as 



r(3r 3 - 2p 2 l) ' 



Kr 3 



3r 3 - 2p 2 l 



(20) 



where K is an arbitrary integration constant. This trajectory is a straight line, and the particle's 
radial velocity is such that it has greater absolute values near the defect. This effect is illustrated 
in Figure 2, where we plot the coordinate r{t) against t (the coordinate system is centered at the 
defect), based on the integration of ([2"U]h 

where a and b are integration constants. We see in Figure 2 the effect of decreasing the absolute 
velocity as particle gets away from the defect (if there was no defect, the curve would be a straight 
line). Notice that as much the particle is away from the defect, its velocity approaches a constant 
value, as one should expect (for r » p, kinematics is the same from a space without defect). 




Figure 2: Radial curve for p = ^P and tq = 50. 



To conclude, let us consider 9 = tt/2 (9 = 0) and <fi ^ 0. The geodesic equations read 

6p 2 l . 2 r(p 2 l + 3r 3 )- 



r + 



r(3r 3 - 2p 2 l) 3r 3 - 2p 2 l 



0. 



(21) 



■■ 4p 2 l- 6r 3 
r(p 2 l + 3r 3 ) 

The last equation can be integrated and we obtain 

Ar A 



rd) = . 



(3r 3 + p 2 l) 2 ' 



(22) 



(23) 



where A is some integration constant. Far away from the defect, we see that <fi is proportional to 



,-2 



Substituting (|23|) into (|21j) . we have 

6p 2 l 



\p 2 l + 3r 3 ) A 2 i 



0. 



r(3r 3 -2p 2 Z) 3r 3 - 2p 2 l (3r 3 + p 2 lf 

The above equation can be numerically integrated, and we can learn that the behavior of radial 
velocity is very similar to the one described by (|20p . which can be seen in Figure 3 (both drawn 
with the help of MAPLE software). 
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Figure 3: Radial curve for p = ^P and tq = 50. 



One can extract an interesting information from (|23|) . together with the natural assumption 
that the absolute velocity decreases in time (as suggested by ([20]) ). As the solution <j> = (A/9)r~ 2 
corresponds to straight lines (geodesies in flat space without defect), the actual derivative <f> given 
by (|23p decays more slowly comparing to the straight line case. This means that the actual path 
must deviate from a straight line, curving to the side of the defect. In other words, the path of a 
test particle is deflected around a defect in a similar way of the gravitational deflection. 



Conclusions 



We study a new kind of defect, which we call ball dislocation, using geometrical methods in linear 
elasticity theory. Whenever the displacement vector (whose discontinuity characterizes the defect) 



is small comparing to natural dimensions of some physical system, the linear elasticity theory is 
suitable, and the formalism of Geometric Theory of Defects can be disconsidered. Moreover, we 
consider a single defect and no other complicated configurations, as a continuous distribution of 
defects (for which the Geometric Theory of Defects is required and well-suited). 

Nevertheless, it is interesting to investigate the formulation of Geometric Theory of Defects for 
our problem. In doing so, we find that a direct (naive) application of this formulation are faced to 
ambiguity problems, in contrast to other kinds of defects (see [6]). The corresponding calculations 
are given in the Appendix. 

Some interesting properties can be seen in the trajectories of free classical particles which follow 
the geodesic equations in the presence of spherical defect. Among the properties of such motion, 
we show that any orbit (around the defect) confined on a sphere is forbidden. The circular orbit 
is a particular case. In the same time we know that circular orbits are permitted in gravitating 
systems; thus, according to at least this feature, the kinematical effects of a defect should not be 
completely identified with gravitational effects. On the other hand, all trajectories are deflected 
near the defect, in an analogous way of gravitating systems. 

One can ask if such a defect could describe some real condensed matter system, where other 
effects than gravity are dominating. In this case, we have a geometric description which mimics 
condensed matter effects from electrodynamics. This question is open, but the present article is 
a first step in studying the issue. It would be natural to identify each atom with a defect, and 
the effects on quantum particles (e.g., Dirac fermions) will be an interesting problem addressed to 
future works. 
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Appendix 

In order to consider the Geometric Theory of Defects, one should start from the induced metric, as 
derived in linear elasticity theory, calculate the corresponding curvature tensors, and identify the 
Einstein tensor with the energy-momentum tensor in the geometric dynamical equations (which is 
the Einstein equations) [3] (see also [B]). 

On can write the line element (1151) in the form 



ds 2 = (1 - vfdr 2 + (r - u) 2 {d0 2 + sin 2 9dcp 2 ) . (24) 

and calculate the components of the curvature tensor, given by 

Rfiupa = d^T upo . - r vp A r AtcrA - (jj, -o- v). (25) 



so 



Rrdr 

ii r( fr r 

R 



-l(r — u) 
-l(r — u) 
-l{r-uf 



5'(r -p) + 
5'{r -p) + 



(1-v) 



(1-v) 



S(r — p) 



5(r - p) 



sin 



-5{r - p) + -77r<5 2 (r - p) 



sm 



_{l-v) K r/ {1-vY 
with 5'(r — p) = d5(r — p)/dr. The components of the Ricci tensor are given by: 

Ru/3 = R^u^p- 

Hence, for the line element (J24D : 

21 



Rfto 



(r — u) 
I 



5'(r-p) 



(1-v) 



5(r - p) 



Rd 



-^{(l-v)(r-u)5'(r-p) + 
(1 — v) A 

+ [v'(r -u) + 2(1 - v) 2 ] S(r - p) + 
+l(l-v)5 2 (r-p)}, 

-3«1 - v)(r - u)5'{r - p) + 

(1 — v) A 

+ [v'(r -u) + 2(1 - vf] 5(r - p) + 
+l(l-v)5 2 {r-p)}sm 2 9. 



The scalar curvature R is given by: 



R = R%. 



Hence 



R 



4/ 



:{(r - u)(l - v)5'(r - p) + 



(r — u) 2 (l — v) 3 

+ [v'(r - u) + (1 - v) 2 ] 5{r - p)} 

2/2 

~7 \2n \2 S2 ( r ~ P)- 

[r — u) z (l — v) z 



(26) 
(27) 
(28) 

(29) 

(30) 



(31) 



(32) 
(33) 



(34) 



Notice that the curvature is non-trivial only in the gluing surface. Moreover, these quantities are 
also ambiguous because of the appearance of the product of 5-function for discontinuous functions, 
and the ambiguity is not cancelled in the calculation of Einstein tensor. 
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